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Abstract 

The two-magnon problem for the frustrated XXZ spin-1/2 Heisenberg Hamiltonian and external 
magnetic fields exceeding the saturation field B s is considered. We show that the problem can 
be exactly mapped onto an effective tight-binding impurity problem. It allows to obtain explicit 
exact expressions for the two-magnon Green's functions for arbitrary dimension and number of 
interactions. We apply this theory to a quasi-one dimensional helimagnet with ferromagnetic 
nearest neighbor J\ < and antiferromagnetic next nearest neighbor J2 > interactions. An 
outstanding feature of the excitation spectrum is the existence of two-magnon bound states. This 
leads to deviations of the saturation field B s from its classical value Bf which coincides with the 
one-magnon instability. For the refined frustration ratio {Jz/Ji] > 0.374661 the minimum of the 
two-magnon spectrum occurs at the boundary of the Brillouin zone. Based on the two-magnon 
approach, we propose general analytic expressions for the saturation field B s , confirming known 
previous results for one-dimensional isotropic systems, but explore also the role of interchain and 
long-ranged intrachain interactions as well as of the exchange anisotropy. 

PACS numbers: 75.30.Ds, 75.30.GW, 75.10.pq, 75.10.Jm 
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I. INTRODUCTION 



We study a spin-1/2 quasi-one dimensional helimagnet with ferromagnetic (J% < 0) near- 
est neighbor and antiferromagnetic (J 2 > 0, IJ2/J1I > 1/4) next-nearest neighbor in-chain 
interactions. In the classical approximation the spins are vectors. In zero magnetic field, 
they form a planar spiral structure (say in the xy plane) with a pitch angle 

cosy? cl = -Ji/4J 2 

between neighboring spins. When a magnetic field is applied along the z-axis, the spin 
moments are inclined toward the z-axis by an angle 

sin0 cl = 8//£J 2 /(4J 2 + J1) 2 , 

where [i = —gfJ>B is the value of the magnetic moment. For fields greater than 

fiBf = (4J 2 + Jx) 2 /8J 2 (1) 

the angle 6 = n/2 and the system becomes "ferromagnetic" (fully polarized uniform state). 

In the quantum case, this high-field ferromagnetic state becomes unstable when the fre- 
quency of a certain excitation mode vanishes. The one-particle instability occurs just at the 
classical field -B^ 1 given by Eq. (PJ. For the collinear antiferromagnet and the obtuse-angle 
helimagnet (if > vr/2) the quantum saturation field coincides with the classical one 1 . In 
contrast, for an acute-angle helimagnet (ip < vr/2) the corresponding saturation field exceeds 
the classical value Bf^il due to the existence of n -magnon bound states below the n-magnon 
continuum (see section II. B of Ref. Q). 

Below we derive an explicit exact expression for the two-magnon Green's function at 
magnetic fields B > B s . It exhibits isolated poles below the two-particle continuum which 
correspond to two-magnon bound states. According to Refs. QJ^, for | J 2 /<A| > |ot c | ~ 0.38 
the two-magnon spectrum minimum determines the saturation field B s > Bf. Our approach 
allows to refine the value of a c , to reproduce their results for B s for the isotropic Ji-J 2 
Heisenberg model and to generalize it to more complex situations of exchange anisotropy 
and interchain interaction as well as of an additional in-chain interaction J 3 . 
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II. THE MODEL AND NOTATIONS 



The Hamiltonian of the model reads 



rrt m <r 



jz nz nz _|_ " r / 0+ q— I c- q+ 



(2) 



where m enumerates the sites in the chain, r determines the nearest (r = ±1) and the 
next-nearest (r = ±2) neighboring sites. We have allowed for an uniaxial anisotropy of the 
exchange interactions. We restrict ourselves to the case of s = 1/2. Then the model given 
by Eq. (J2J) can be applied to undoped edge-shared chain cuprates^. Here the spin operators 
may be expressed via the hard-core boson operators b m 
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(3) 
(4) 
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where the square brackets stand for the commutator, and m denotes the site index. The 
ferromagnetic state corresponds to the vacuum state b \ FM) — b |0) = 0. Then the Hamil- 
tonian (0) can be rewritten as 



H = Hq + H,, 



inti 



<m ^ rrt r- 



(6) 
(7) 



^ mt = 2? J ^ m "'" 



2 

r i 



+r- 



(8) 



The transverse part of H © defines the one-particle hoppings in H Q (J7J), the Ising part 
contributes the interaction (jHJ) and on-site energy value u . 

We shall study the one- and two-particle excitation spectra of the Hamiltonian given 
by Eq. (JHJ) which will be obtained from the singularities (poles and branch cuts) of the 
corresponding retarded Green's functions (GF) 



Gi ja (k,uj) 



Ak,i\At a 



(9) 
(10) 
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where 



{(X\Y)) = -i dte 1 ^-^ 
Jt' 

1 



X(t),Y(t>) 



A 



k.l 



k/2+qOk/2-q 



1 Y,e- lk{m+l/2) b m b m+l . 

i V rn 



(11) 



The expectation value (...) denotes the ground state average, the time dependence of 
an operator X(t) is given by X(t) = e ltH Xe~ ltH , and the Fourier transform of b m reads 
b q = iV -1 / 2 I] m exp(— zgm)6 m . N denotes the total number of sites. 



III. THE ONE-MAGNON SPECTRUM 



The equation of motion for the hard-core boson operators (jHJ) reads 

1 _ 



4 K 



b m ., H 



Uobm + 



m I u m+r J-r^m+rbm 



(12) 



For the ferromagnetic ground state, the terms proportional to n do not contribute to the 
one-magnon GF (jOJ). This means that the usually infinite hierarchy of equations of motion 
including higher order Green's function is cutted exactly and closed rigorous expressions for 
all n-magnon Green's function can be obtained in principle. In particular, the one-magnon 
GF becomes simply 



GV>(q,u,) = ((bM 



UJ — U! 



sw 



where 



sw 



1 



^o+oY, J r V ex P inr) 



(13) 



(14) 



* r 

is the free spin-wave dispersion. 

The dispersion has a minimum at the helical wave vector q = ip cl /a, where a = 1 
is the lattice constant. 

For the anisotropic Hamiltonian (J2J) it is convenient to define 



a = J% y /J* y , Aj EE J Z JJ* V . 



Then 



cos<£> cl = — l/4a, 



(15) 



(16) 
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and for magnetic fields values smaller than 



fiBf = J{ + jf + r 2 y + 



+ A 2 



-, , (4a+l) 2 



uiq™ becomes negative. Evidently, in the isotropic case Ai 



(17) 



A 2 = 1, Eq. JIZD reduces to 



Eq. ©. 



IV. THE TWO-MAGNON GREEN'S FUNCTION 

The operator A^i (|11|) annihilates a pair of particles, separated by the distance / and 
moving with total quasimomentum k. A two-particle bound state manifests itself by an 
isolated pole of the two-magnon GF (TMGF) (110(1 . The hard-core condition (jSJ) demands 
A k}0 = 0. We see also from Eq. (fTTj) that A k j = A k -i ■ The time evolution of A k j, I > is 
given by the relation 



d j 
at 



A kjl , H\ = 2u A Kl + 
1 



<N 



-ik(m+l/2) 



£ 



r i <"m+r u m 



Jm+l 



+ A' r 



(18) 



Using the commutation relations (@J), and the symmetry J_ r = J r we rewrite the operator 
A k i in the normal form 



A'm,l ~ E 
r 

rrfim+l 



+ 



Again we note that for the ferromagnetic ground state, the terms containing the operators 
n do not contribute to the GF and as discussed in the previous section corresponding higher 
order GF vanish exactly. Then, within the subspace of two-particle excitations above the 
ferromagnetic ground state, we may write rigorously 



A k ./, H 



(2u + Jf) A Kl + (1 - S h0 ) £ J r V cos —A Kl+r . 



(20) 



Thus, the problem of calculation of the TMGF (|10|) is equivalent to the impurity problem 
for the one- dimensional tight-binding-like Hamiltonian 



H tb {k) = f + V, 



(21) 



f = 2loq \ m ) ( m \ + '^2\m + r) t r (m\ , 

m m,r 

V = \ m ') £ m' (m'\ , 

III' 

where 

t r (k) = Jf v cos — , m' = 0,r, e = 00, e r = J*. (22) 

Let us note that the infinite value of £0 is the result of the hard-core constraint given by Eq. 
(J5J). The periodic part T results from H (JJ|), and H int (JHJ) defines the changes of the on-site 
energies on impurity sites. 
It is easy to see that 

G lta (k,u) = (<f )l \(u-H tb y 1 \<p a ), (23) 

where \^) = + \-j)) /y/2, j = I, a. 

In a standard way, we will use the identity 

(u-H tb y l = (uj-ty' (24) 
+ (uj-ty'v {u-H tb y l 

for the solution of the impurity problem in the real space. After some algebra we obtain 

G ltl (k, u) = l — , (25) 

~ Jl + r (o) j zr (0) 

n (0) . ^1, 2 j 2 Lt 2,1 
^ + 1 r (P) jz 

where Gfl is the GF of non-interacting hard-core bosons (H = Hq) 





^9l9a /nr\ 

gi+a + gi-a , (26) 

go 



1 cos ql 

N ^ , , (, ,sw 1 , ,sw v ! ' ' ' 



J k/2+q ^ ^k/2-q 

V. THE TWO-MAGNON BOUND STATES AND THE SATURATION FIELD 

In the derivation of the exact expression for the Green's function ([25 )1 - ([27 )1 we have 
used the mathematical equivalence of the Heisenberg model (J2J) at high fields with the ID 
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impurity problem (|21|). But the obtained explicit expressions have a rather complicated 
form. Fortunately, the physics of the same ID impurity problem helps also considerably in 
its further analysis. 

The branch cut of Gi^{k,uj) (J25j) is defined by the continuous part of the spectrum of 
Htb(k) given by Eq. (|21j) that corresponds to the two particle continuum of the starting 
Hamiltonian @. Its boundaries may be found from the spectrum of the periodic part T 

E(k, q) = 2[u + h{k) cosg + t 2 (k) cos2g] = u s k f 2+q + u%j%_ q . 

Since ti(k) < for all k, we have 

E{k,qi) < E{k,q) < E{k,ir), \k\ < fa, (28) 
E(k,0) < E(k,q) < E{k,n), fa < \k\ < fa, (29) 
E(k,0) < E(k,q) < E(k, q x ), fa < \k\ < tt (30) 

where 

cosgi = -ti(£0/4t 2 (A;) 
E(k, gi ) = 2uo -tl/4t 2 -2t 2 , 



y/128a 2 + 1 + 1 . . . , 

fa = 2arccos < tt/2, tAfa) At 2 {fa) = -1, 

16|a| 

fa = 2arccos — 28 ^, + , 1 — - > n/2, ti(fa) / At 2 {fa) = 1. 
lo|a| 

As we will see below, the point k = n has a special meaning. For this value of k the 
nearest-neighbor hopping £1(71") = 0, and the Hamiltonian T describes two non-interacting 
linear chains (i.e. the sites with odd or even numbers m) with a hopping t 2 {^) = —J 2 V inside 
each chain. The account of the hard-core constraint (JHJ), i.e. the introduction of Vo = |0)e (0| 
with Eq = 00 entering the Hamiltonian H t b{k) given by Eq. (|2*T|) , does not influence the 
chain with odd sites, but the chain with even sites m is broken into two independent semi- 
infinite chains. Now, it is easy to account for the rest of terms in the impurity Hamiltonian 
V , because E\ = J\ and e 2 = J 2 affect different chains. The GF G 2i2 (tt,uj) has a particular 
simple form. It can be obtained e.g. by the recursion method 

G 2 , 2 (it,uj)= L-(2u + J^)-(J^) 2 G^(n,uj)r 1 = -^- ^ — y , (31) 

J J 2 [z — ^2 — T {z)\ 
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where the dimensionless energy 

z = (u- 2u )/J% y (32) 

is introduced, and r(z) = (z — V z 2 — 4^ /2 = 1/ [1 — r(^)] is the local Green's function on 
the first site of the unperturbed semi-infinite chain in dimensionless units. A simple analysis 
of the expression (j31)l shows that besides the branch cut in the interval —2 < z < 2 of 
the real axis, G^^ 71 "?^) may have an isolated pole. The pole exists for A 2 > 1 above the 
continuum, i.e. a bound state exists for the easy-axis anisotropy of the next-nearest neighbor 
exchange. The expression for G\ j i{ , k,u) is more complicated than (j2,2(tt, ^) 

G ltl (ir,u) = [lo- (2uj + J{ - J* y ) - J?t(z)\- x (33) 

= Jf [z - Ai/a + 1 - t(z)) ■ (34) 
Note that for any acute angle helimagnet (a < 0) a bound state should exist below the 
continuum. Indeed, the condition Gi\(k,Ub) = gives 



For fe 7^ vr, a bound state exists, too^. In Refj2j, the isotropic version (Ai = A 2 = 1) of 
the Hamiltonian (J2J was considered. There A. Chubukov has found that the dispersion of 
the two-magnon bound state Zb(k) exhibits a minimum at k = n for \a\ > \a c \ ~ 0.38. The 
latter number will be refined below. 

Based on extensive numerical work^ for finite chains, and qualitative discussion in Ref. 
0, we strongly believe that the two-particle bound state defined by (jBEJ) is the excitation 
with the lowest energy per flipped spin in the system for this parameter regime. The 
absolute dominance of two-magnon states manifests itself by AS* 2 = 2 steps of the calculated 
magnetization curves M(H) at high fields. Only for a <■ 0.4 steps with AS Z = 3 are observed 
(see Fig. 1 of Ref. Iff). But we admit that from a formal point of view, the full rigorous solution 
should also include the analysis of the problem for the arbitrarily n-magnon bound-states 
(n > 3) in a similar way as done here. However, the corresponding calculations are rather 
cumbersome and particular examples (n = 3, 4) are left for future consideration. 

Then the quantum saturation field is determined by the condition that the two-magnon 
energy vanishes (i.e. the two-magnon instability of the field-induced ferromagnetic ground 
state). This way, the central result of the present work yields 

u b {ir, B s ) = 2/1(5 - B s ) = J^z b (rr) + 2cu = 0, (36) 



or 



fj,B s 

Which gives explicitly 



B„ 



E^-^(vr). 



2 J?{ Jf + J% y ) - ( Jf f - 2 J| Jf 



(36') 



2/i(Jf 



2// 



2 "Jf) 

2(A 2 + l)-(A 1 /a) 2 -2A 1 /a 



(37) 



1 - Ai/a 

For the isotropic case, this result was first obtained in Ref. |2J by solving an integral equation 
which results from a summation of a sequence of ladder diagrams. From our straightforward 
derivation it is clear that the result is exact within the adopted two-magnon approach, as 
also pointed out in Ref. 0. 

In order to find the parameter range where the value of Zb(Tc) (see Eq. ([35))) yields the 
minimum of the bound state dispersion z b {k) , we consider the expressions (|25l V f|27|) in the 
vicinity of k = ir, z = z b (k). After straightforward calculations given in the Appendix we 
obtain 



z b (k) « z b (ir) + 



(k 



7T 



2m 



eff 



where 



2m e //(a, Ai, A 2 



1 Aj_ 

2 4aAf 



a a 
+ a 



2Ai 



(2a- Ai)A 2 



2(a-A 1 y 4A X (a - A x ) [a (1 + A 2 ) - AJ 



(38) 



(39) 



The dependence of the right-hand side of Eq. ()39|) on the inverse frustration ratio | J\j J 2 | 
is shown in Fig. ^ for different values of the nearest-neighbor exchange anisotropy Ai . 
The dependence on the second neighbor exchange anisotropy is weak in the vicinity of the 
isotropic point A 2 = 1. We see that m e ff is positive for large values of frustration J 2 ^> | J\\ 
and changes the sign at some critical value a c , where the dispersion minimum is transformed 
to a local maximum. For A 2 = 1 the condition l/2m e //(a c , Aj, 1) = reduces to a cubic 
equation for a c and we have 



a c (Ai) = Ai 



2 J R(Ai)cos 



0(Ai) + 2vr 



(40) 



where 



i?(Ai) = 



arccos 



54R 3 
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0.5 1 1.5 2 2.5 3 3.5 4 
Inverse frustration ratio 1 /oc=| J-, |/J 2 

FIG. 1: The inverse effective mass (|39|) as a function of frustration ratio and anisotropy Ai = 
Jf j J* y for the J\ — J2 model. 

The critical frustration dependence on anisotropy is shown on Figure El • In the isotropic 
case, we have R(l) = —2/3, 0(1) = arccos (7/16), and 

a c (l) = 1{ 2 cos + 1^ « -0.37466105983527, (41) 

which refines |a c (l)| ~ 0.38 calculated before^. If Ai — > 0.5, then a c from Eq. (j4Tj|) diverges, 
i.e. the effective mass m e ff(Ai = 0.5) becomes negative for all frustration values. 

The Figure EH shows the quantum and classical saturation field dependencies on the pa- 
rameter values of the ID isotropic J\ — J2 model. We have chosen J2 as the unit of energy. 
We see that the quantum effect is most pronounced for frustration values \J%/ Ji\ ~ 1- 

In the region 0.25 < \a\ < \a c \ the minimum Zb(k) for the isotropic model shifts into 



the points k = 2qo = 2 arccos(— l/4a). The authors of Ref. argue that in this case 
the saturation field is determined by bound states of three and/or more magnons-; such a 
situation is out of the scope of this paper. For the anisotropic J\ — J2 model, there is a 
third possible scenario. The one- and two-particle instabilities occurs at different /c-points. 
Then, it is possible to have the minimum of z^k) at k = rr , which is higher in energy than 
the lowest boundary of the continuum (f28j) z c (2q ) = z c (2q ). This happens e.g. for the 
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FIG. 3: The two-particle {ijlB s ^/J2 - solid line) and one-particle (/j,Bf 2 /J2 - dashed line) values of 
the saturation field as a function of frustration ratio | J1I/J2 for the isotropic J\ — J2 model. 



11 



easy-plane nearest-neighbor anisotropy values 



1 + VI + 16a 2 . 
Ai < Ai a = v —— , A 2 = 1 

8|a| 

In Fig. 0] the dependence of the saturation fields on the anisotropy of the J\ ex- 
change is shown. The former is important for edge-shared cuprates 7 . We see that for 
Ai = (1 + v / 17)/8 » 0.64039 the lines £ s (Ai) and 5 s cl (Ai) do intersect. At the same time 
l/2m e ff(l, 0.64039, 1) ~ 0.16178 > 0. The unexpected at first glance result that the clas- 
sical curve apparently reaches, then overwhelms the quantum result can be explained by 
the reduced attractive ferromagnetic interaction due to the anisotropy, i.e. a weakening of 
the two-magnon "glue". Below this value of Ai the saturation field is determined by the 
one-particle instability, like in the XY model (Ai = A 2 = 0). Thus, the intersection is not 
related to a strange quantum versus classical behaviour, but to the competition between 
one- and two-particle instabilities. 

It is interesting that a strong easy-axis anisotropy can diminish the saturation field, and 
at the point 

Ai, = (-a) 



1 + y/3 + 2A 2 

the field B s vanishes, i.e. the system's ground state becomes ferromagnetic^. Note also the 



region 



1 + A 2 + 



a 2 



< A x < Ai o 



where Bf = 0, B s ^ 0. Here the classical fully polarized state is destroyed by quantum 
fluctuations. A possible ground state for the system in this parameter regime may be a 
collinear state with period 4 described in Ref. yl or a dimer nematic state, predicted in Ref. 
for the isotropic model in a high field. 



VI. ADDITIONAL INTERACTIONS 

An advantage of our approach is the possibility to apply it to more complex situations 
which occur naturally when real chain compounds are considered. 

Indeed, it is easy to realize that the exact mapping of the two-magnon problem onto the 
effective tight-binding Hamiltonian (|21|) is not restricted to ID and to the J\ — J 2 model. 
We may generalize the Hamiltonian H © including into summation over r more distant 
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C\J 




FIG. 4: The quantum {^B s ^/ J2 - solid line) and classical (fiBf 2 /J2 - dashed line) saturation fields 
for the anisotropic J\ — J 2 model, J\ = J 2 = A 2 = 1. 

neighbors in chain direction. This will introduce additional impurities and hoppings in 
the effective Hamiltonian (J21|) . Moreover, we may consider also 2D or 3D systems. Then, 
the site indices m as well as r in ((7J), (jSJ) becomes vectors with corresponding changes in 
the effective Hamiltonian 1)21)1 . It is straightforward to obtain the TMGF (123)) . but the 
expression becomes cumbersome. Here we will apply our general approach to answer the 
question how is the quantum effect for the saturation field in the J\ — J 2 model modified by 
some additional interactions often present in real compounds. 

First, we include a small third neighbor in-chain interaction J 3 . Such a term may appear 
as a result of the spin-phonon interaction in the antiadiabatic regime, when the exchange 
constants Jj <C hw p h the characteristic phonon frequencies, and the spin-phonon interaction 
is strong*- It is expected to be small J3 <C | J\ \ , J 2 and antiferromagnetio^. Below, the 
subscript 2(3) refers to the J\ — J 2 and the Ji — J 2 — J 3 model respectively. In this section, 
for the sake of simplicity we give only formulae for the isotropic case J = J z = J xy . The 
minimum of the one-magnon spectrum ()14j) gives the value of the helicoidal wave vector and 
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the classical value of the saturation field 

-J2 



cos g ,3 



,cl 



Jf — 3Jz( J\ — 3J3 



6J3 

A 3J 3 
cosg ,2 ^1 - -y- 

^ 2 + ^fi-- 



1 



/i \ 4a 16a 1 -' 

we recall that a < 

The two-magnon Green's function (|33|) has the form 



(42) 
(43) 
(44) 



Gi,i(7r,o;) 



a; — (2co> + Ji — -h 



J 2 

■J 2 



~1 -1 



^-(2 Wo + J 3 ) + J 2T (^; 



(45) 



For small J3 values, the bound state energy and the saturation field varies linearly with J3 

a(2 — a) 



B 



^b,3 ~ ^6,2 — ^3 
J 3 



+ 1 



a 



s,3 



2/i 



1 + 



a(2 — a) 
(1-aY 



(46) 



where the values Ub,2, and B Si 2 are given by Eq. (J37|) . The slope of dependence on J 3 is 
larger than for -B Sj 3. It means, that positive J3 suppress the quantum effect. The difference 
of saturation fields in quantum and classical cases becomes smaller. 

As the simplest example for a two-dimensional system we consider a 2D set of chains 
parallel to the x-axis coupled in perpendicular direction with the strength J±. Then the 
one-magnon dispersion becomes two-dimensional 



^Jd = ^g x W iD + J± (cos q y a - 1) . 



,sw 



(47) 



From this expression one readily obtains 



^2D = VBlD + ^-L+ I J-L 



(48) 



i.e. in this approximation the ferromagnetic interchain interaction does not affect the satu- 
ration field, whereas in the antiferromagnetic case it is enhanced by 2J±. In general, such a 
correction is especially important near the quantum critical point for ferromagnet-helimagnet 
transition a m (4 + 9 J3/ J2}~ 1 , where the ID saturation field by definition vanishes. Eq. ()48jl 
should be understood as a lower bound for the saturation field near the critical point. The 
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account of quantum fluctuations will lead to slightly higher values of B s according to Ref. 



For an arbitrary k-point, the GF (J23j) is found from the solution of a system of three 
linear equations, but along the line k = (7r/a,k y ) the system reduces to a single equation 
which gives 



where the two-dimensional spectrum (jTfj) should be used in the expression for the non- 
interacting GF (J26|) . The dispersion of the isolated pole and the two-particle continuum 
boundary are shown in Figs. 03 and El For small J± (Figure EJ) one observes a well separated 
bound state. Here, the absolute minimum of the continuum occurs at k = (2g ,0) and at 
k = (0, 0) its energy z c = — 2.45 Ji exceeds the minimum of bound state dispersion given 
by Z(,(0,7r) = — 2.61576J2- For strongly coupled chains such with J± = \ Ji\ = J 2 , the pole 
position becomes very shallow (Figure EJ) and it becomes clear that such a local minimum 
exceeds the minimum given by two independent (one-magnon) excitations (z c = — 4.25J 2 
for the parameter set shown in the caption of Fig. EJ). In the general 3D problem, one may 
expect that even the bound state itself may disappear. 

At variance with the classical case given by Eq. (}4*8|) . the solution of Eq. ()49|) yields for 



i.e. in this case the saturation field is sensitive to both sign of interchain interaction. With 
the increase of J± one finally reaches a critical value, where z^ = z c and beyond the "one- 
magnon" derived Eq. (|48|) should be used instead of a "two-magnon" one like Eq. (|5()j) . 
Thus, the quantum effects are maximally pronounced in the ID case, just as the localization 
for the equivalent impurity problem. 

More complex interchain interactions derived from band structure calculations and inelas- 
tic neutron scattering datap^ and an application to chain cuprates will be given elsewhere. 

VII. CONCLUSION 

We have shown that the internal motion of a two-magnon pair on a ferromagnetic back- 
ground is equivalent to the motion of a single particle described by an effective tight-binding 
Hamiltonian. This Hamiltonian is not translationally invariant. It models the hard-core 
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(49) 



J± |< J2 



fiB s = liB Si1D + J ± + 0{J 2 JJ 2 ) 



(50) 
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FIG. 5: The two-magnon bound state energy Zb{-K,k y a) = (w&(7r , k y a) — 2uq)/J2 (solid line) and 
the boundary of the two-magnon continuum z c (tt, k y a) (dashed line) as a function of the quasi- 
momentum value in the y direction for |Ji| = J2, the interchain interaction is chosen as J± = 
O.IJ2, J2 being the unit of energy. For these parameters, the absolute minimum of continuum is 
z c (2q ,k y a) = -2.45J 2 > z b 

boson constraint (jSJ) by an infinite on-site energy at the site with zero coordinates, and 
each exchange J r introduces the on-site energy e r and the hopping term t r = J r coskr/2. 
Remarkably, this mapping procedure can be applied to problems at arbitrary dimension. 
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k 

y 

FIG. 6: The same as in previous figure for | J\\ = J2 = J±. The absolute minimum of continuum 

is z c (2qo,k y a) = — 4.25J2 < z\> 

The two-magnon Green's function is found exactly by analogy with the impurity problem. 
The two-magnon excitation spectrum is found from poles and branch cuts of the Green's 
function. For the quasi-one- dimensional helimagnet with ferromagnetic nearest neighbor 
and antiferromagnetic next-nearest neighbor interactions a bound state of magnons exists. 
This leads to deviation of the quantum saturation field B s from the classical value. 

The derived expression for the saturation field B s (exact within the two-magnon ap- 
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proach) provides a constraint for competing exchange interactions. Such a constraint may be 
useful in fitting thermodynamic properties such as the magnetic susceptibility x(T) and the 
magnetic specific heat c p (T). In general, high-field magnetization measurements M(H,T) 
yield an important information concerning the exchange integrals in novel materials. Com- 
bined with the analysis of other experimental data, this knowledge may be very helpful to 
elucidate the relevant microscopic model for an acute-angle helimagnetic system (i.e. having 
a pitch angle < 7r/2 at zero magnetic field). 

In this work we studied the lowest energy of excited states, i.e. the position of the isolated 
TMGF poles. The obtained Green's function given by Eq. (110(1 contains the information 
about the whole spectrum that is necessary for the calculations of physical properties for 
concrete materials. Various application to edge-shared compounds will be considered else- 
where. 
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APPENDIX A 

Here we give details in derivating the Eqs. (JH8|) and (J2HJ). 

The energy of an isolated pole of GF given by Eq. is the root of the equation 

~ J l + r (0) T z r (0) = °' ( Al ) 

(0) , ( - T l,2 J 2 tj 2,l 



1 r (0) jz 
1 — { ~ T 2,2 J 2 



It depends on the k value via the dependence of the hopping parameters t r in T ()2Tjl . We 
denote k = n — k, and expand the left-hand-side of Eq. (jAljl up to terms oc k 2 . Then 

h = Jt y sin £ « Jf^, t 2 = -Jl y cos k « -J? (l - y ) • 
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Note that the GF given by Eq. flUJ) 

9i{k, u) = g-i(k, uj) = (R\ (u - ff 1 \R + I) , 
obeys the equation of motion 

(w - 2lu ) g l = <^ + tx (gi-i + g t+1 ) + t 2 (g^ 2 + gi+2) ■ 



(A2) 



We will calculate go(k, u),gx(k, uj) directly from Eq. (J27|) and use Eq. (|A2J) for the calculation 
of the other g\ involved into Eq. (jAljl . We begin with 

dQ 



(A3) 



uj — 2 (ujq + tx cos Q + t 2 cos 2Q) 

The denominator of the integrand is nonzero for uj outside the spectrum of T. After the 
substitution r = tan(Q/2) a straightforward calculations give 

1 



£0 



8 J 2 y cos kJ (4g + p 2 ) (q + p — 1) 
p-2 - V4g + p 2 



p- 2 + V4g + p 2 



q + 1 — v / 4gTp 2 " 



g + 1 + V4g + p 2 . 



where 



P 



2sin(K/2) 



2 — 2 COS K 



w — 2uj 
J* y 



4a cos k 4cos/t o 2 

Expanding this expression around the point (k — ir, z = Zb(ir)), we obtain 



1 



go 



Analogously we obtain 
gi(k,uj) = — 



2 A lA2 4A 2 a 2 - 3a 2 + 3Aia - Aj 
1 K " (J _ r 2A?(2«-A 1 ) 



(A4) 



cos Qc/Q 



27r J-7T — 2 (cuq + tx cos Q + £2 cos 2Q) 



2J 2 y ^/(4q + p 2 ) (q + p-1) 
ko (a — Ax) 2 



2q + p- y/Aq + p 2 2g + p + v^g + p 2 



g + 1 - V4g + p 2 V<? + 1 + V4g + p 2 



2Jf A?(2a-Ai)" 
Now, using Eq. (jA2|) . we have 



(A5) 



g2 + go 



1 ~ yf+i 
2 J?" 



1 - «■ 



2a A 2 (2a 2 - 4A x a + A 2 ) - 3 (a - A] 
4Af (2a - AT) 



\3l 



(A6) 
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Substituting the above expressions (|A4[) . ()A5|) . ()A6|) into Eq. (j2Hj). we obtain 



^1,1 



-i 



2Jf 



1 - / z ~ 2 

V z+2 



1 + K 2 



a (2a 2 - 4Ai« + A?' 
A?(2a-A0 



4Af (2a - Ai 



and 



^(o) _ r (o) 

Gg(7T,^(7r)) 



k a 

2jf A 2 ' 

a 



J? (Ai - a) ' 

In the neighborhood of the point k = tt, z = z^tt) Eq. (|A1|) may be rewritten as 



-4 + Gf3 



0. 



(A7) 



(A8) 
(A9) 



1- JiG$(ir,z b (n)) 

The substitution of Eqs. (|A7|) . (|A8jl and (|A9|) into this equation allows to solve it with 
respect to z and to obtain finally Eq. (J38j) . 
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